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Abstract 



The hadronic parameters of pseudoscalar (r/ c ) and scalar (xc) charmonium are determined at 
finite temperature from Hilbert moment QCD sum rules. These parameters are the hadron mass, 
leptonic decay constant, total width, and continuum threshold (so). Results for so(T) in both 
channels indicate that so(T) starts approximately constant, and then it decreases monotonically 
with increasing T until it reaches the QCD threshold, sth = 4toq, at a critical temperature 
T = T c ~ 180 MeV interpreted as the deconfinement temperature. The other hadronic param- 
eters behave qualitatively similarly to those of the J/ip, as determined in this same framework. 
The hadron mass is essentially constant, the total width is initially independent of T, and after 
T/T c ~ 0.80 it begins to increase with increasing T up to T/T c ~ 0.90 (0.95) for \c (%)> and sub- 
sequently it decreases sharply up to T ~ 0.94 (0.99) T c , for Xc (Vc), beyond which the sum rules 
are no longer valid. The decay constant of Xc at first remains basically flat up to T ~ 0.80 T c , 
then it starts to decrease up to T ~ 0.90 T c , and finally it increases sharply with increasing 
T. In the case of r] c the decay constant does not change up to T ~ 0.80 T c where it begins a 
gentle increase up to T ~ 0.95 T c beyond which it increases dramatically with increasing T. This 
behaviour contrasts with that of light-light and heavy- light quark systems, and it suggests the 
survival of the rj c and the Xc states beyond the critical temperature, as already found for the J/ip 
from similar QCD sum rules. These conclusions are very stable against changes in the critical 
temperature in the wide range T c = 180 — 260 MeV. 
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1 INTRODUCTION 



The extension of the QCD sum rule method [T] to finite temperature was first discussed in [2], 
and subsequently developed successfully to obtain the T-dependence of hadronic parameters, and 
study the phase transitions of chiral-symmetry restoration and quark gluon deconfinement [3] - 
[7j . Further evidence for the validity of this extension to finite temperature was given in [8] . It has 
been established in this framework that hadronic spectral functions experience substantial modi- 
fications with increasing T. In particular, in the case of light-light and heavy-light quark systems, 
their leptonic decay constants, weak and strong couplings, hadronic widths, and form factors de- 
velop a prominent T-dependence consistent with deconfinement at some critical temperature T c . 
The advocacy of a substantial temperature dependence of hadronic masses has been refuted by 
results from these applications. In fact, these masses hardly change with T. If sq is the squared 
energy beyond which the hadronic spectral function merges into a continuum, well represented 
by perturbative QCD (PQCD), it has been found that s (T)/s (0) ~ (0|gg|0)| T /(0|gg|0}| T =o [3J. 
The vanishing of the quark condensate signals chiral symmetry restoration, while the vanishing 
of the continuum threshold sq(T) is an indication of quark-gluon deconfinement, as resonances 
disappear and the smooth spectral function is determined purely by PQCD. The gluon conden- 
sate plays essentially no role here due to its peculiar behaviour, i.e. it hardly changes until T 
becomes very close to T c . However, this will not be the case for heavy-heavy quark systems due 
to the absence of the light-quark condensate. Hence, a qualitatively different temperature be- 
haviour of charmonium is to be expected. In fact, we have shown recently [7] that the J /tp total 
width initially increases with increasing T up to some temperature T* < T c after which it begins 
to decrease beyond T c , while the leptonic coupling remains approximately constant up to T*, 
and subsequently it increases with temperature also beyond T c . This behaviour provides strong 
evidence for the survival of charmonium in the vector channel beyond the critical temperature 
T c , in agreement with lattice QCD results [9]. The T dependence of sq(T) in the J/ip channel 
is qualitatively similar to that for light-light and heavy-light quark systems, i.e. sq(T) decreases 
monotonically with increasing T until it reaches the QCD threshold Sth = ^ m Qi with tuq being 
the quark mass. Beyond this point there is no longer a solution to the FESR as the integration 
range vanishes. The dynamical origin of this behaviour lies with the gluon condensate as well 
as with the T-dependence of the PQCD vector correlator in the space-like region, the so called 
current-quark scattering term. In fact, there is a distinct interplay between these two terms with 
increasing temperature. 

At this stage it must be pointed out that in the framework of QCD sum rules the critical tem- 
perature for deconfinement, referred to above, is just a phenomenological parameter. It is the 
temperature at which the resonance couplings and the continuum threshold approach zero, and 
the widths increase sharply, for light-light and heavy-light quark correlators. Hence, it need not 
coincide numerically with e.g. the critical temperature obtained in lattice QCD, which is defined 
differently. In fact, results from QCD sum rules lead to values of T c somewhat different from 
those from lattice QCD. Hence, comparisons between different frameworks should be made in 
terms of the dependence of parameters on the ratio T/T c . 

In this paper we extend the analysis of [7j to charmonium in the scalar (xc) and the pseu- 
doscalar (r/ c ) channels. We find a temperature behaviour qualitatively similar to that of J/ip, 
i.e. a monotonically decreasing so(T), suggesting deconfinement at a critical temperature, a total 
width which initially is independent of T, and then it increases with T up to T/T c ~ 0.90 — 0.95, 
beyond which it begins to decrease dramatically, and a leptonic coupling which starts constant, 
then begins to decrease for Xc (or already starts to increase for the r] c ), and finally shoots up 
sharply in both cases. The hadronic masses remain essentially constant. In these two channels 
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there is no contribution to the correlator in the space-like region (scattering term), so that the 
gluon condensate is responsible for all of the thermal behaviour of these hadronic states. Due to 
this circumstance, the critical temperature for deconfinement, T c , corresponds to that at which 
the gluon condensate vanishes. Given the uncertainty in the value of T c we consider the wide 
range T c ~ 180 — 260 MeV. Qualitatively, the results remain very stable in this range. Slight 
quantitative variations, however, do not affect the conclusions. Since sq(T) approaches the QCD 
threshold sth = 4?tiq as T approaches T c , the integration range in the FESR approaches zero, 
and the method is no longer valid. Nevertheless, the dramatic decrease of the width, and the 
increase of the decay constant near T c suggests the survival of the r\ c and \ c beyond the critical 
temperature for deconfinement, in agreement with lattice QCD results [9]. 



2 QCD SUM RULES 

We consider the correlator of the heavy-heavy quark (pseudo)scalar current at finite temperature 

n (5) (g 2 ,T) =i j d^xe^ 9(x ) « |[J( 5) (x) , J ( f 5) (0)]| » , (1) 

where J( 5 )(x) =: Q{x)TQ{x) :, T = 75 or Y = 1 for the (pseudo) scalar case, and Q(x) is the heavy 
quark field. The vacuum to vacuum matrix element above is the Gibbs average 

<< A ■ B »= exp(-E n /T) < n\A ■ B\n > /Tr(exp(-H/T)) , (2) 

n 

where |n > is any complete set of eigenstates of the (QCD) Hamiltonian. We shall adopt the 
quark-gluon basis, as this allows for the standard QCD sum rule program at T = to be seamlessly 
extended to T 7^ [4]. The QCD sum rules to be used are the Hilbert moments 



1 / d \ i N + l ) 1 f°° ds 

^(QlT) = n (5)( Q3,r)| Q2 ^ = - I {s + Ql)(N+2) ImT K s,T) , 

(3) 

where N = 1,2, and Q 2 Q > is an external four-momentum squared to be considered as a free 
parameter |10j . Using Cauchy's theorem in the complex s-plane, which is equivalent to invoking 
quark-hadron duality, the Hilbert moments become Finite Energy QCD sum rules (FESR), i.e. 

<Pn(QI,T)\had = <Pn{Qo, t )\qcd , (4) 

where 

I r s o(T) j s 

Vn(Q ,T)\had = - ( s + Q2)(jv+ 2 ) ImU(s,T)\ HAD , (5) 
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1 f so{ - r > ds 

<Pn(Qo,T)\q C d = ~ , n 2\( N +D I™>RpQCd(s,T) + <pn(Qo,T)\np ■ (6) 

The imaginary part of the (pseudo) scalar correlator in perturbative QCD (PQCD) at finite tem- 
perature, Im II(5)((/ 2 , T), involves two pieces, one in the time-like region (q 2 > 4toq), Im (q 2 , T), 

which survives at T=0, and one in the space-like region (q 2 < 0), Im H^(q 2 ,T), which vanishes 
at T=0. A straightforward calculation of the pseudoscalar correlator in the time-like region, at 
T ^ 0, and to leading order in PQCD, gives 



-Iml4 a) (^) = -W v(uj 2 ) 



TT 



8tt 2 



2 **(£) 



4m 2 Q ) , 



(7) 



where [«(ur)] 2 = 1 — Auiq/uj , mg is the heavy quark mass, q = u — q = oj in the rest frame 
of the thermal bath, and np(z) = (1 + e z )~ 1 is the Fermi thermal function. The quark mass is 
assumed independent of T, which is a good approximation for temperatures below 200-250 MeV 
The equivalent expression for the scalar correlator is 



llmU^(oj,T) = u 2 [ v (oj 2 )Y 



TT 



1 -2n F 



\2T 



6(uj 2 - 4m 2 ,) 



(8) 



Both spectral functions are strongly exponentially suppressed. In the so called scattering region 
{q 2 < 0) the scalar QCD spectral function is given by 



-ImU (s \q 2 ,T) 

TT 



16vH 



np 



|q|x + uj 
2T 



np 



q|x — ui 
2T 



dx 



q 

2T 



dz 



Integrating by parts, the above equation becomes 



0) 



-ImU (s \q 2 ,T) 

TT 



2tt 2 Q \ I-' 3 



0J 

|q| 



Q ( m Q\ . 

~r 2 np {-T) +2 



znp{z) dz 



(10) 



where lim| q |^ |q| v = 2 uiq has been used above. Invoking the limit 



UJ 



lim - — — = - 8(uj 



|q|->° q' 



(11) 



the QCD spectral function, Eq.(10), vanishes since lini| q |^ q 2 = lim q |^ \ (uj 2 — |q| 2 ) 5(uj 2 ) = 

| uj 2 5(uj 2 ) = 0. For the pseudoscalar correlator the QCD spectral function in the scattering region 
is 
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rip 



2T 



rip 



\q]x — to 
2T 



dx 



12) 



In the rest frame of the thermal bath this contribution vanishes, i.e. 



(13) 
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Figure 1: The ratio so(T)/sq{0) as a function of T/T c , for r/ c (curve (a)), and \c (curve (b)) 
corresponding to T c = 180 MeV. For T c = 260 MeV the break-point is at T/T c ~ 0.55, instead of 
T/T c ~ 0.80. 



This leaves the gluon condensate as the leading temperature dependent term in the OPE. The 
moments for this term have been calculated at T = in |10j . After taking into account that our 
definition, Eq.(3), differs from that in [10J by one additional derivative, the results for the scalar 
correlator are 
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<Pn(Q 2 ,T)\ np 



27 2( 7V+1 )AH 1 (iV + l)(iV + 2)(iV + 3)(iV + 4) 

"8^ 2 (4m2 (2iV + 5)H (27V + 7) 



F ( iV + 2,--,iV + -;p ; 3(iv + 4) 



2 / 19 

" WiV + 2,-,iV + -;p 



$,(14) 



where F(a,b,c,z) is the hypergeometric function, , £ = p = j^g, and 



$ = 



4vr 2 1 



9 (4m 



2 \2 

Q> 



a 



—G 2 



TT 



(15) 



where {{^fG 2 )) stands for the temperature dependent gluon condensate. The result for the 
pseudoscalar is 



<P5n(Q 2 ,T)\ np 



3 2( Ar+1 )iV! 1 (N + 1)(N + 2)(N + 3)(iV + 4) 



8vr 2 (Am 2 Q )( N + 1 ) (1 + f) 



N+2 



(2iV + 3)H (2iV + 4) 



F[N + 2-~,N+ 7 -;p 



- — - — -F ( N + 2.--.N+ -;p 
(N + 4) V 2 2 



$(16) 



The thermal behaviour of the gluon condensate was obtained in [12J; to a good approximation 
this result can be written as 



—G 2 



— G 5 

TT 



9(T* - T) + 



1 - X. 



1 - ^ 



9(T - T* 



(17) 



where T* ~ 150 MeV is the breakpoint temperature where the condensate begins to decrease 
appreciably, and Tc ~ 180 MeV is the temperature at which ((^■C 2 ))^ = 0. It should 
be mentioned that this lattice QCD determination has rather large errors, so that this value 
of T c is affected by a large uncertainty. We have changed the value of T c in the wide range 
T c ~ 180 — 260 MeV, and found no qualitative change in the results, as discussed in the next 
section. This simple parametrization of the lattice data can be improved using a smooth func- 
tion without changing the conclusions to be reached with Eq.(17). At T=0 the gluon conden- 
sate has been extracted from data on e + e~ annihilation, and r-lepton hadronic decays |13j : 
i^fG 2 ) = (0.05 ± 0.02)GeV 4 . We do not consider a gluonic twist-two term in the OPE [T4] as it 
is negligible, as in the case of vector (J/VO charmonium [7]. In fact, after calculating the Wilson 
coefficients corresponding to the (pseudo) scalar correlator, the contribution from this term is 
roughly one order of magnitude smaller than the scalar gluon condensate at T = T c . 



Turning to the hadronic representation of the (pseudo) scalar correlator, we follow the standard 
procedure and parametrize it in terms of the ground state resonance, i.e. the r\ c (xc), followed 
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Figure 2: The ratio of the r? c mass M(T)/M(0) as a function of T/T c for T c = 180 MeV. A 
similar result is found for the Xc as well as for T c = 260 MeV. 



by a continuum given by PQCD after a threshold sq > M 2 , where M is the hadron mass. This 
ansatz is an even better approximation at finite temperature, as sq(T) decreases monotonically 
with increasing T in all systems analyzed so far. Considering first the zero width approximation, 
the hadronic spectral function is given by 



- ImU(s,T)\ H AD = -ImIl(s,T)\ R Es 0(s - s) + -ImIl(s,T)\pQCD 0(s - s ) 

7T TT TT 

= f(T)M 2 (T)6(s-M 2 (T)) +-ImU(s,T) a e(s-s ) , (18) 

TT 



where s = q 2 = ui 2 — |q| , and the leptonic decay constant is defined as 



<0|J (5) (0)|ff(fc) >=fM 2 , 



(19) 
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Figure 3: The ratio of the n c width r(T)/r(0) as a function of T/T c corresponding to T c = 
180 MeV. For T c = 260 MeV the break-point is at T/T c ~ 0.55, instead of T/T c ~ 0.80, so that 
the curve becomes roughly twice as broad. 



where H(k) stands for rj c {xcj- Next, considering a finite (total) width the following replacement 
will be understood 



S(s - M 2 (T)) =^> const - — , 0/ - , (20) 

V V >' {s-M 2 (T)) 2 + M 2 (T)T 2 {T) ' 1 ; 

where the constant is fixed by requiring equality of areas, e.g. if the integration is in the interval 
(0 — oo) then const = M(T)T(T)/ir. To complete the hadronic parametrization one needs to con- 
sider the hadronic scattering term due to the current scattering off heavy-light quark pseudoscalar 
mesons (D-mesons) in the thermal bath. In the case of the rj c , the pseudoscalar current couples 
to an odd number of pseudoscalars, so that the scattering term is effectively a negligible higher 
order (one-loop) effect [3]. For the scalar case it is easy to show that the hadronic scattering 
spectral function vanishes in the rest frame of the thermal bath, as this term is not singular as 
|q| — s- 0, and u — >■ 0. 
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Figure 4: The ratio of the Xc width r(T)/r(0) as a function of T/T c corresponding to T c = 
180 MeV. For T c = 260 MeV the break-point is at T/T c ~ 0.55, instead of T/T c ~ 0.80, so that 
the curve becomes roughly twice as broad.. 



3 RESULTS 



We start the analysis at T = 0, and use the moments, Eq.(3), to reproduce the experimental 
values of the mass for an input width of the r] c and Xc [15j . The leptonic decay constants are not 
known experimentally, so that they will be predicted from the sum rules, having in mind that they 
are expected to be a few hundred MeV. Ultimately, their precise value is not important, as what 
matters is the behaviour of their ratio /(T)//(0). Beginning with the zero-width approximation 
it follows from Eq.(5) that 



MQ 2 )\had _ MQ 2 )\had , , 

MQ 2 )\had ^{Q 2 )\had ' 1 j 

which turns out to hold with extremely good accuracy given the small widths of the r\ c and the 
Xc- Using Eq.(4) this leads to 



8 





1.5 




1.4 




1.3 




1 .z 


✓ — V 




o 


1.1 


P 


1.0 








0.9 




0.8 




0.7 




0.6 




0.0 0.1 



Figure 5: The ratio of the leptonic decay constant / as a function of T/T c for the r\ c (curve (a)), 
and the Xc (curve (b)) corresponding to T c = 180 MeV. For T c = 260 MeV the break-point is at 
T/T c ~ 0.55, instead of T/T c ~ 0.80. 



Vi(.Qq)\qCD = ¥2{QI)\qCD 
¥2(QI)\qCD V2,{QI)\qCD ' 



(22) 



which depends only on the QCD information, involves the two unknowns so and Qq, and provides 
the first equation to determine this pair of parameters. It is very important to stress at this point 
that for consistency sq must be determined from the sum rules themselves, rather than importing 
it from some other analysis. The second equation can be Eq.(4) with e.g. N = 1. Following 
this procedure, and starting at T = we find for r\ c that so(0) = 9.1 GeV 2 , Q\ = 0, and 
r(0) = T\ EX p = 27 ± 3 MeV give M(0) = 2.9 GeV, to be compared with the experimental 
value M(0)\ EX p = 2.98 GeV. For the Xc we find that s (0) = 11.6 GeV 2 , T = T\ EX p = 
10.2 ± 0.7 MeV, and Ql = yield M (0) = 3.2 GeV, to be compared with the experimental value 
M(0)\exp = 3.41 GeV. The leptonic decay constants are predicted to be / = 183 MeV for rj c , 
and / = 190 MeV for % c , in agreement with expectations. The above results are fairly insensitive 
to changes in and N in the range = — 10 GeV 2 , and N = 1 — 6 (for larger values of N 
there is no solution for sq(0)). A posteriori, and at finite temperature, varying Qq and N in that 
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wide range leads to no appreciable changes in the results and conclusions. 

Next, at finite T Eq.(21) continues to hold with extreme accuracy even if the width were to grow 
substantially. In fact, for a tenfold increase in the width, Eq.(21) remains valid within a fraction 
of 1%. Solving for sq(T) the results are shown in Fig.l for T c = 180 MeV, where curve (a) is 
for the pseudoscalar channel, and curve (b) for the scalar channel. The end point of each curve 
corresponds to so(T c ) = 4:mQ ~ 6.8 GeV 2 , i.e. the QCD threshold value which then makes the 
FESR integrals vanish. The break-point at which sq(T) starts to decrease occurs at T/T c ~ 0.80, 
for T c = 180 MeV (T/T c ~ 0.55 for T c = 260 MeV). Since in the present case there is only one T - 
dependent driving term, it is not possible to extend the analysis beyond the critical temperature. 
For the J/ip there are two driving terms, the PQCD scattering term proportional to T 2 , and 
the gluon condensate defining a critical temperature [7]. This feature allowed the analysis to 
extend beyond T = T c . Once again, sq(T) must be obtained from the QCD moments themselves, 
rather than using some temperature parametrization from some other application. Otherwise, 
results could be seriously inconsistent. The result for the mass of the r\ c is shown in Fig. 2; a 
similar behaviour is found for the Xc- This result confirms once again the view that the hadronic 
mass is not a meaningful order parameter, as it hardly changes with temperature; in some cases 
decreasing slightly, and in other cases increasing slightly with increasing T. Conceptually, there 
is no particular reason for the mass, being the position of the spectral function pole on the real 
axis, to be a signal for deconfinement. In contrast, the total width and the coupling do provide 
such a deconfinement signal, if the width increases and the coupling decreases with increasing 
temperature. In view of this result for the hadronic masses, we keep them constant in the sequel 
to simplify the analysis, and discuss later the consequences of allowing a slight T-dependence. 
Results for the widths are shown in Fig. 3 for r/ c , and in Fig. 4 for Xc- Qualitatively, they resemble 
the result for J/ip. Quantitatively, the rj c width is essentially independent of the temperature up 
to T/T c ~ 0.80 (T/T c ~ 0.55 for T c = 260 MeV), after which it increases substantially by a factor 
nine until T/T c ~ 0.95 (T/T c ~ 0.87 for T c = 260 MeV) where it decreases dramatically by more 
than a factor two at T/T c ~ 0.99 (T/T c ~ 0.97 for T c = 260 MeV). Beyond this temperature 
there is no longer a solution to the FESR. The behaviour of the Xc width is qualitatively similar 
to that of the r\ c . Quantitatively, though, the Xc width increases by almost a factor eighteen 
until a slightly lower T/T c ~ 0.90 (T/T c ~ 0.75 for T c = 260 MeV), and then it decreases by 
a factor two at the end point T/T c ~ 0.94 (T/T c ~ 0.85 for T c = 260 MeV). The behaviour 
of the leptonic decay constants is shown in Fig. 5 for rj c (curve (a)), and for Xc (curve (b)). For 
the rjc, the temperature at which the coupling starts to increase dramatically is close to the 
temperature at which the width begins to decrease substantially, i.e. T/T c ~ 0.95 (T/T c ~ 0.97 
for T c = 260 MeV). The end point is the same as for the width, i.e. T/T c ~ 0.99 (T/T c ~ 0.97 
for T c = 260 MeV). A similar correlation is found for the Xc, i-e. the coupling reverses behaviour 
at T/T c ~ 0.90 (T/T c ~ 0.75 for T c = 260 MeV), and the end point occurs at a similar 
temperature as for the width, i.e T/T c ~ 0.94 (T/T c ~ 0.85 for T c = 260 MeV). The origin 
of this correlation between the width and the coupling, and the fact that there is no longer a 
solution for the width and coupling (end point) below T c is easy to understand from Eqs.(4)-(6). 
For instance, in the zero-width approximation the square of the coupling is proportional to the 
right hand side of Eq.(6), i.e. to the sum of the PQCD integral and the non-perturbative (gluon 
condensate) moment, the latter being negative (see Eqs. (14) and (16)). At some T < T c the 
monotonically decreasing PQCD integral becomes equal to the non-perturbative moment so that 
beyond this temperature f 2 < 0, and there is no longer a real solution. In the case of the widths 
the end point corresponds to the temperature beyond which they become negative. We comment 
in closing on the validity of the simplifying approximation M(T) ~ M(0). Allowing for the 
hadronic mass to change (decrease) near T c produces a slight reduction, together with a shift to 
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the right, of the peak in Figs. 3 and 4, but preserving the overall qualitative behaviour of the 
widths. Something similar happens with the couplings, i.e. the temperature at which they start 
to increase substantially shifts slightly toward T c . 



4 Conclusions 

We have determined the temperature dependence of the hadronic parameters of charmonium 
in the (pseudo) scalar channels, ij c and Xc using Hilbert moment QCD sum rules. The results 
confirm the expectation that heavy-heavy quark systems behave differently from light-light and 
heavy-light quark hadrons, due to the absence of the light quark condensate in the OPE. Since 
in the space-like region there is no contribution to the PQCD nor to the hadronic (pseudo) scalar 
correlator, this leaves the gluon condensate as the leading term in the OPE, and thus solely re- 
sponsible for the temperature behaviour of rj c and Xc- The continuum threshold sq(T) was found 
to follow closely the T-dependence of the gluon condensate. The leptonic decay constant and the 
width behave initially as in light and heavy-light hadrons, i.e. they both are initially constant, 
and after certain temperature the former decreases and the latter increases with increasing T. 
However, beyond a certain temperature this behaviour reverses, with the decay constant increas- 
ing and the width decreasing with T up to some temperature just below T c beyond which the 
sum rules are no longer valid. In the (pseudo) scalar channel it is not possible to carry on the 
analysis beyond T c , as in the case of J/ip, because the gluon condensate is the sole driver of the T 
- dependence (there are no scattering terms as in the vector channel). Nevertheless, the behaviour 
of the decay constant and the width of r] c and Xc near T c suggests their survival beyond T c , as 
found for the J/ip in this same framework [7], and in agreement with lattice QCD results [9]. 
Experimentally, a possible way of confirming the survival of charmonium states beyond T c would 
be to observe them in heavy ion collisions at high beyond a certain threshold value. 
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